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Abstract 

o : 

CN ■ This is a technical work about how to evaluate loop integrals appearing in one loop 

nonplanar (NP) diagrams in noncommutative (NC) field theory. The conventional wisdom 

Ph ■ says that, barring the ultraviolet /infrared (UV/IR) mixing problem, NP diagrams whose 

planar counterparts are UV divergent are rendered finite by NC phases that couple the 
loop momentum to the external ones p through an NC momentum p M = Q^ v p v . We show 
■ that this is generally not the case. We find that subtleties arise already in the simpler 

case of Euclidean spacetime. The situation is even worse in Minkowski spacetime due to 
its indefinite metric. We compare different prescriptions that may be used to evaluate 

© ; l oop integrals in ordinary theory. They are equivalent in the sense that they always 

yield identical results. However, in NC theory there is no a priori reason that these 

^ ■ prescriptions, except for the defining one that is built in the Feynman propagator, are 

Q_i' physically justified even when they seem mathematically meaningful. Employing them 

can lead to ambiguous results, which are also different from those obtained according to 
> ■ the defining prescription. For p 2 > 0, the NC phase can worsen the UV property of loop 

^ ■ integrals instead of always improving it in high dimensions. We explain how this surprising 

phenomenon comes about from the indefinite metric. This lends a strong support to the 
point of view that the naive approach is not well-founded when time does not commute 
with space. For p 2 < 0, the NC phase improves the UV property and softens the quadratic 
UV divergence in ordinary theory to a bounded but indefinite UV oscillation. We employ 
a cut-off method to quantify the new UV non-regular terms. For p 2 > 0, these terms 
are generally complex and thus also harm unitarity in addition to those found previously. 
As the new terms for both cases are not available in the Lagrangian and in addition can 
be non-Hermitian when time does not commute with space, our result casts doubts on 
previous demonstrations of one loop renormalizability based exclusively upon analysis of 
planar diagrams, especially in theories with quadratic divergences. 

PACS: ll.10.Nx, ll.10.Gh 
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1 Introduction 



Quantum field theory on noncommutative (NC) spacetime may be formulated on ordi- 
nary commutative spacetime through the Weyl-Moyal correspondence, with the point- wise 
product of field operators in the action replaced by the Moyal star product, 



where x, y are ordinary commutative coordinates and 9^ v is the NC parameter matrix 
characterizing the NC spacetime, = i9^ u . 

The star product induces a momentum-dependent phase in vertices of Feynman di- 
agrams in momentum space, which could potentially improve the convergence property 
of loop diagrams. However, there always remain a group of diagrams, i.e., planar ones, 
that are only modified by a global phase that does not touch a loop momentum. For 
these diagrams, the loop integrals are exactly the same as their ordinary counterparts 
and thus the ultraviolet divergences persist [1, 2]. More interesting are the diagrams that 
do couple external momenta to loop momenta via the NC phases. These nonplanar (NP) 
diagrams, as usually called, exhibit some exotic features that are unfamiliar in ordinary 
field theories. For instance, they break unitarity of the S matrix when time is involved in 
noncommutativity [3]. Nevertheless, there seems to be no doubt that they are ultraviolet 
finite at one loop level, barring the ultraviolet /infrared (UV/IR) mixing problem that 
happens for exceptional kinematic configurations [4]. Then, at least at one loop level, 
the fate of renormalizability of an NC field theory rests on the consistent removal of UV 
divergences in planar diagrams. Based on this understanding, it has been extensively 
demonstrated in the literature that one loop renormalizability is guaranteed for the NC 
version of popular commutative theories; see Refs. [5, 6], for instance, for the first few 
examples about U(l) and U(N) gauge theories. 

In this work, we examine this statement about finiteness of one loop NP diagrams. 
The issue turns out to be much more delicate than expected. An subtlety appears already 
in NC theories on Euclidean spacetime: whether an NP integral converges or not can 
depend on how we evaluate it. As this essentially mathematical subtlety does not seem to 
be solvable by mathematics, we suggest in section 2 a physical motivation to fix it in the 
context of quantum field theory. The situation on Minkowski spacetime is more involved 
due to its indefinite metric. To prepare for later discussions, we review in section 3 some 
prescriptions that may be used to evaluate loop integrals in ordinary theory. They are 
equivalent in the sense that they always yield the same answer. Then, we present our 
main results in section 4, based on the defining prescription used in Feynman propagator. 
These results are quite different concerning convergence property from the conventional 
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ones. The other prescriptions are examined in section 5. We show that they are generally 
not equivalent to the defining prescription any longer, even if they sometimes seem to 
make sense mathematically. For an illustration of physical implications of the results 
in section 4, we compute in section 6 the one loop NP contributions to the self-energy 
and vertex in real y? 4 theory and discuss briefly some new points about unitarity and 
renormalization. Finally, we summarize our results in the last section. 



We consider in this section the simpler case of Euclidean spacetime where the integrand 
of a loop integral is real and positive definite before the NC phase is introduced. The 
integral of an NP diagram at one loop in n dimensions is essentially of the type, 



where m 2 > 0, k 2 = k^k^ > 0, p 2 = p^p^ > and k • p = k^p^. Here p^ = 0^ u p u with p 
a linear combination of the external momenta. For a non-tadpole diagram, one first uses 
Feynman parameters to combine propagators and makes a shift in the loop momentum. 
In that case, m 2 is actually a linear combination of the internal masses squared and the 
quadrature of the external momenta, and those Feynman parameters have to be integrated 
over upon finishing the above k integral. The relevant points for us are that m 2 is still 
positive definite and that introducing Feynman parameter integrals does not modify UV 
properties of the original loop integral. For a convergent integral it is allowed to shift 
the loop momentum. For a nonconvergent one, we assume that it has been regularized. 
For integrals that involve the phase e lk ' p , we assume that it is legitimate to ignore the 
sin(A; • p) term odd in k. In the strict sense this is valid only for convergent integrals. For 
non-convergent ones and their regularized versions, we will take their principal values in 

/+A r+oo 
dx f(x) = for an odd function f(x) even if / dx f(x) does 
-A J-oo 

not exist. Finally, the rotational symmetry in p is employed in the above notation. 

We only need to evaluate A™ since A$ is related to Ag^ by a derivative in m 2 . This 
is obvious for convergent integrals. For non-convergent ones, it is also true as long as the 
regularization method used to define them is independent of m 2 . This will be the case 
in our cut-off method. We will not discuss the UV/IR mixing problem, thus we assume 
p 2 ^ in our discussion. Then, with rescaling \fp 1 k ii = x^, we have, 



2 Euclidean spacetime 
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with p being the unit vector of p. To illustrate the subtleties hidden in the original integral 
A™, we employ several methods to work out E n (a). Then, we discuss how to manipulate 
these integrals in the context of field theory. 



2.1 Mathematical subtlety in E n (a) 
2.1.1 Method 1 

Without loss of generality, we choose the axis Xi to be in the direction p. We finish first 
the x\ integral to obtain, 



cxp 



a 



E»{a) = n Uidxj) J ==L= J , (4) 



/= 2 1 v/E^f + a 2 
which is obviously convergent. Finishing the angular and radial integration, we have, 



E n (a) = -Qn^a 71 ' 2 r dx x^ ( 
2 Jo 



V^ + T ( 5 ) 
= (27r)5at- 1 ir 1 _ f (a), 

where fig is the solid angle in £ dimensions and K„(z) is the modified Bessel function of 
order v [7], so that 

A™(p 2 ,m 2 ) = (27r)-t(m 2 )t- 1 (^rfp^j 2 K X -^ yJmP^j (6) 

According to this method, the above result holds in any dimensions. This coincides with 
the conventional one in the literature. 



2.1.2 Method 2 

In this method, we reverse the integration order in method 1 to finish first the subintegral 
over components perpendicular to the direction p. It is obvious that for n > 3 the 
subintegral diverges so that the integral itself is not well defined, and that for n = 2 
the result in method 1 is reproduced. These two methods make explicit reference to the 
direction defined by p. 

2.1.3 Method 3 

Noting the rotational symmetry of the integral E n (a), we finish first the angular inte- 
gration in n dimensions. Although any method of angular integration must lead to the 
same answer, a judicious choice of angles facilitates this. We choose X\ = jjcos^.! with 
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x = ^Jx^x^ and Q n -\ G [0, 7r] so that, 



/o x A + cr 

where J„(z) is the Bessel function of the first kind of order v. The leading term of J v (x) 
in the large x > limit is, 

J v (x) — > y — (~ sin a; or ~ cos a?) (8) 

V TTX 

Thus the above integral converges only for | — | — 2 < 0, i.e., for n < 5, in which case 
we obtain, 

E n (a) = (2rr)^a^- 1 K^ 1 (a) (9) 

the same result as in method 1 upon using K-„(z) = K u (z). For n > 5 however, the 
integral is not well-defined in this third method. 

2.1.4 Method 4 

The last method is the one using Schwinger parametrization that is popular in the liter- 
ature. In Euclidean space, the denominator of the integrand in E™ is positive definite, 
hence it is justified to exponentiate it as follows, 

E n (a) = J d n x e ix -P da e~ a{x2+a2 \ (10) 

where the cosine factor is replaced by the phase in the sense of principal value mentioned 
earlier if convergence is not guaranteed. 

To make the exponentiation useful for evaluating the integral, we have to interchange 
the integration order of x and a, which however is justified only for (absolutely) conver- 
gent, or regularized if not, integrals. Suppose we take this for granted to proceed further. 
Then, complete the quadrature in k, make a complex shift x — > x+ that is permitted 
because the integrand is analytic in each complex plane of x. Finishing the x integral 
yields, 

E n (a) = Trf I" da a^e~ aa2 -^ = (2tt)? a^K^ia), (11) 
jo 2 

which holds seemingly for any dimensions as in method 1. 
2.2 Discussions 

The above calculations tell us that NP integrals in two dimensional Euclidean spacetime 
are indeed unambiguously finite, barring the UV/IR mixing problem. However, in three 
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and higher dimensions different computational methods yield different results for the same 
integral. In this respect we do not count method 4 as a justified one, because the very 
use of Schwinger parametrization presumes that the integral is well-defined and thus it 
cannot say anything when the integral is not well-defined. On the other hand, the other 
three methods are completely justified but still yield different results. This looks quite 
puzzling. This is a new feature of NP diagrams in NC field theory that does not appear 
in ordinary field theory and that seems to have not been realized previously. 

The answer to the dilemma is simple: the integral A™ as originally defined in eq. (2) 
is not mathematically well-defined in three and higher dimensions. To define it prop- 
erly, we should also specify its integration path or order. Method 1 amounts to preferred 
integration over the component in the p direction. As the phase is active only in this 
direction, this manipulation earns the best convergence property for the remaining com- 
ponent integrals. Method 2 is just the opposite. It starts with the components that are 
independent of the phase. The subintegral thus corresponds to an ordinary integral in 
(n — 1) dimensional Euclidean spacetime and divergence starts to appear at this stage 
for n — 1 > 2. This method gives the worst convergence property. Method 3 starts with 
integrating spherically over directions in n dimensions for a fixed radius and integrating 
over the radius later. This produces a mild convergence property. 

The result of method 1 can also be understood from another angle, which explains 
clearly why the similar subtlety does not appear in ordinary theory. Take E 5 (a) as an 
example which converges according to method 1 but does not according to methods 2 and 
3. We employ the NC phase cosrr! to do integration by parts repeatedly, reducing the 
UV divergence degree by one each time, until we arrive at an integral that converges by 
power counting. Since the denominator is positive definite in Euclidean spacetime, the 
surface terms vanish. For E 5 (a), we need to do so three times, 



E b (a) = 24 / n( rf;r i) / dx i sinx i 

m J 



xi 2x\ 



(EL^ + « 2 ) 3 (£L^ + a 2 ) 4 



(12) 



Since the above is now convergent, we can interchange the integration order to finish Xj^\ 
first and x\ afterwards, yielding the same result as obtained in method 1. 

In the context of field theory, we integrate over momentum fluctuations of virtual 
particles in loops. When the above subtlety appears, mathematically speaking, there is 
no a priori reason as to which way of integration we should take. But physically, method 
3 is much better motivated than methods 1 and 2 for the following reason. A single 
Feynman diagram usually involves several NP terms that share the same denominator 
(i.e., m 2 ) from propagators but differ in the NC phases (i.e., p M ) from the decomposition 
of the product of vertices. In method 3, we treat all terms in the same diagram (actually 
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also all diagrams) on the same footing with no preference to any of the p^s. In methods 
1 and 2 however, we accumulate contributions of fluctuations in a manner that varies 
according to their preferred directions specified by p M 's, which seems rather unnatural. 

Let us consider the critical case of n = 5 according to method 3. To regularize the 
integral, we impose a spherical cut-off on the original integral, k^k^ < A 2 , which amounts 
to a cut-off x^x^ < k 2 on the integral E 5 (a) with k = h-yfp 1 . We obtain 

E 5 (a) = 87r 2 / dx^- — -(^-^-cosrr) 

Jo x 2 + a z \ . x J 

= —8tt 2 sin k + (2ir)^a^K^[a) H , 



(13) 



(14) 



where the dots stand for terms vanishing in the limit A — > oo. Thus, 
Al(p 2 ,m 2 ) = -(47r 3 )- 1 (p 2 )-Isin(^A^ 

+ (2tt)-3 (m 2 )i 2 Kz (yfrrtp) , 

which has a new term compared to the conventional result. It oscillates but is bounded 
in the UV. In even higher dimensions, the new terms are no longer bounded though 
oscillating in the UV. For instance, we quote the result for n = 6 without giving details, 



A\(p 2 ,m 2 ) = -2(27r)~lv / A(p 2 )-^cos(^Ayp 2 -^7r 
+ (27r)- 3 (m 2 ) 2 (v^V)~ 2 K 2 (v^V) 



(15) 



From the above derivation, we observe that Ag converges for n < 45 + 1, oscillates at 
n — AS + 1, and diverges oscillatorily for n > 45 + 1. Although convergence of NP loop 
integrals is improved over their commutative counterparts, new terms will appear in high 
enough dimensions that are non-regular in the UV and cannot be removed by sending the 
cutoff A to infinity. 

In any cut-off method, there is always a question concerning when to impose the cut-off 
as the method generally breaks rotational or Lorentz invariance. In the above calculation, 
the cut-off is imposed after the original integral has been cast into the form of eq. (2). 
This guarantees rotational invariance with respect to external momenta contained in m 2 
and p 2 . The cut-off method is still not unique, of course. Using the spherical cut-off as 
suggested by the physical motivation, any manipulation, including method 2, will yield 
the same answer, though method 3 is the easiest in practice. Mathematically speaking 
and without affecting rotational invariance in external momenta, however, we are not 
forced to employ the spherical one. For instance, we may impose two independent cut- 
offs, Ami, on the loop momentum components perpendicular or parallel to p. The above 
mathematical subtlety then appears in the form that the convergence property depends 
on the order of taking the limits of Aj_ — > oo and An — > oo. 
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3 Minkowski spacetime: review of ordinary theory 



Minkowski spacetime has an indefinite metric that produces nontrivial analytic properties 
in Green's functions. The same metric will cause special difficulties in NC theory. Since 
the Feynman propagator will necessarily be defined on the complex plane, an NC phase 
on the real axis ceases to be a pure phase any longer. This will modify the analytic 
properties of Green's functions significantly. Without taking care of this, we could make 
manipulations, especially analytic continuation, that cannot be justified. This indeed 
happens, as we will show in the next sections, when we take for granted that some 
prescriptions equivalent in ordinary theory are equally justified in NC theory. To better 
see the matter, we review in this section these manipulations for ordinary loop integrals. 

3.1 An apparent problem 

The one loop integrals in ordinary theory in n dimensional Minkowski spacetime can be 
cast into the form, 

f d n k 1 

iB^m 2 -ie) = J Qir)" (k* - m* + ie) s ' ™ " 2 ' 5 ~ 6 = ° + ' 

(16) 

with the help of Feynman parameters and free shift of the loop momentum. Now m 2 , as 
a linear form of the internal masses squared and the quadrature of the external momenta, 
can be positive (below threshold), zero (on threshold) and negative (above threshold). For 
n > 25, Bg is divergent in the UV region of k where \k 2 \ is also large, i.e., with k being 
deep space-like or time-like. This divergence is familiar to us. It is essentially Euclidean 
in the sense that it can be figured out by power counting as in Euclidean space without 
considering complications from the indefinite metric. It is handled by regularization. 
However, there seems to be another divergence in the above integral, from the UV region 
where \ko\ and \k\ are equally large but \k 2 \ is small, i.e., the deep light-like region of k. 
It is possible only for an indefinite metric and can be seen more clearly by changing the 
integration variables to x± = k ± \k\ upon finishing the angular integration, 

■Tin, 2 • \ fi n-lo f°° A f°°A (x+ + X-) n ~ 2 

(27r) n J-00 Jo (— 4x + x- — m 2 + ie)° 

where the right-hand side is not well-defined for any 5 in the region where one of \x±\ 
is large but is small. This apparent divergence is often overlooked justifiably by 

practitioners of loop calculations because, as detailed below, it does not actually appear in 
ordinary theory. However, it is questionable that it can also be ignored beyond ordinary 
field theory, and particularly in NC theory where time does not commute with space. 
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3.2 The defining prescription 

It is instructive to see how the apparent problem is avoided in ordinary theory. For 
this, we trace back to the definition of Feynman propagator that is the only ingredient 
appearing in the above integral. 

iA F (x-y) = (OITXxMj/)) |0) 

= 9(x - y ){0\<p(x)<p(y)\0) + % - x o ){0\<p(y)<p(x)\0) (18) 
= 6(x - y )iA(+\x -y) + 9(y - x )iA^ + \y - x), 

where <p is a neutral, free scalar field with mass m and |0) is the free vacuum state. It has 
been expressed in terms of the iA^ function which itself is defined as the commutator 
between the positive- and negative-frequency parts of the field. In four dimensions, it is 

iA^(x) = / ^e~ ikx , k = \Jk 2 + m 2 > m (19) 

[2tiY J 2ko 

Using an integral expression for the step function, the above can be cast into a four- 
dimensional integral, 

iAp(x) = I d ' k A - 1 e~ lk - x (20) 
V ; Jc F (2tt) 4 k 2 - m 2 K ' 

where Cp is a contour in th e comple x k plane. It starts from — oo o n the rea l axis, moving 
below the pole at k = —\Jk 2 + m 2 and above the pole at k = +\Jk 2 + m 2 , to reach +oo 
on the real axis, then moves along the semi-circle of infinite radius in the upper or lower 
half plane according to xo < or x > respectively, to return back to — oo on the real 
axis. The above is usually written in a concise form using Feynman's +ie prescription, 

iAp(x) = /|L F _i r --e- (2!) 

where e is an infinitesimal real positive number and is the source of the +ie appearing 
in the loop integral Bg. Although the contour notation is suppressed in the above, it 
should be understood as being enforced. The two poles have been displaced properly, 
so that the contour now extends over the whole real axis of ko from — oo to +oo and 
returns to — oo via the semi-circle described above. Namely, it is understood that the ko 
integral is always meant to be finished first by contour integration. This is necessary to 
keep the above expression consistent with the definition of Ap(x): positive-frequencies 
propagate only for positive time x$ > while negative-frequencies propagate only for 
negative time Xq < 0. The latter cannot be guaranteed without enforcing the contour 
integration. Furthermore, if the contour integration were not taken into account properly, 
Ap and thus B% would be afflicted with the apparent problem discussed in subsection 
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3.1. This point is often not emphasized in the commutative theory since the prescription 
is equivalent to the one of Wick rotation that is easier to implement. 

With the above understanding of Feynman propagator in its manifestly covariant form, 
our loop integral becomes using residue theorem, 



iBx[m z - ie) = / - — r — - / =; 

SK ' J (27T)"- 1 Jc F 2tt [jfeg - (k\+ m 2 - ie)} 5 



i( 1)' r( ^~ 1} f ^ (V I m 2 ie)hs 



(22) 



where indeed there is no apparent divergence problem that we worried about earlier. Here 
the contour C F can be closed either in the upper or in the lower half plane. Finishing the 
angular integration in (n — 1) dimensions, then \k\, and using product theorem for the T 
function, we obtain, 



n—3 



,r(25-i) a n ^ r°o x - 



m 2 -ie)z- d ^ ^ (23) 



2^[r(5)] 2 (27r)-i v ; r ( 5 "l) 



r (s - -) 
= (-i)* y m 2j wH^- i ey 

where we have assumed n < 25 to make the integral convergent. Divergent integrals can 
be handled as usual. 



3.3 Wick rotation 

As long as we know the apparent problem in subsection 3.1 does not actually appear, we 
can employ other prescriptions that can lead to the same result for loop integrals. The 
only thing that we should take care of is that we must not change the analytic properties of 
the original integral. Amongst those prescriptions Wick rotation is most appealing since 
it can be understood as analytic continuation from Minkowski to Euclidean spacetime. 
Namely, the ko integral can be considered as being integrated over the imaginary axis 
in the complex plane of ko- Since the poles of ko are located in the second and fourth 
quadrant, the rotation has to be counter-clockwise, ko — > e l ^k n , with k n to be integrated 
over the real axis from — oo to oo. The integral can be readily finished since we are now 
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in Euclidean loop momentum space, 

iB^-ie) = e^-^)[^L 

5 v ; J (2?r)" (k 2 E + m 2 - ief 

= i -l V — / cfc- -7 (24) 

identical to the result according to the defining prescription. This prescription is indeed 
easier to implement than the first. 

3.4 e parametrization 

The third prescription is built in the Schwinger parametrization based on the positive 
infinitesimal e, here named e parametrization for simplicity. Relying on e > 0, it is 
legitimate to exponentiate the denominator in B$, 

"*<"*'-*> = ml^r^'^^ (25) 

To employ the exponentiation, we have to interchange the integration order. However, this 
is permitted only for a well-defined integral. Here we are not concerned with Euclidean 
divergences by power counting, since they can always be handled by regularization. In- 
stead, we worry about the potential problem discussed in subsection 3.1. Interchanging 
the integration order presumes that the integral is indeed free of such a danger, and thus 
has to be justified. This we will do later. 

The k integral can now be finished trivially to yield, 

B r i(m 2 -it) = (An)-^i- s e-^ n ^ J™ daa 6 ~^ l e- a{t+tm2) (26) 

Note that the Euclidean UV property of the original integral has been translated into the 
property of the parameter integral at a — > 0. The apparent problem that was avoided 
in the defining prescription by finishing first the k contour integration or in the Wick 
rotation by integrating over the imaginary axis of k , is now hidden in the behaviour at 
a — > oo which is tamed by e > 0. The a integral can be finished straightforwardly to give 
the same answer as the previous prescriptions. But this is not yet what we want, because 
the above interchange of order has to be justified. In addition, it is rather disturbing that 
convergence at a — > oo relies on an infinitesimal e. We recall that e was introduced only 
as a short-cut for routing around the poles and its magnitude should be irrelevant. In 
fact, these two matters are the two faces of the same coin. Namely, the feasibility of the 
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e parametrization relies on the possibility that dependence of convergence on e can be 
relaxed. 

The integrand of the a integral is analytic in its complex plane, with a cut along the 
negative real axis when n is odd. For m 2 > 0, we can rotate the integration path to the 
negative imaginary axis, a — > e~ 1 ^ f3 with f3 G [0, oo), because the exponent has a negative 
real part in the fourth quadrant to offer an exponential decay in the contribution from 
the quarter-circle at infinity: 

B%{m 2 -ie) = (ATt)-^(-l) 6 -^— H d(3 (3 s -^- l e-^ m2 ~ le) (27) 

V / 

whose outcome is the previous result in the case m 2 > 0. For m 2 < 0, we rotate the 
integration path to the positive imaginary axis, a — > e 1 ^ (3. Again, the exponent provides 
a sufficient decay for the contribution of the quarter-circle in the first quadrant. The 
result is, 

B%(m 2 -ie) = (Any^e-^^ J™ dp pt-^e-^™ 2 ^ (28) 

which indeed reproduces the previous result for m 2 < by noting that (— \m 2 \ — ie)%~ 5 = 
e l7T ( s - -) \m 2 \%~ s . The case m 2 = is not a problem either. Remember that we have to 
integrate over Feynman parameters hidden in m 2 . It is dangerous only when m 2 = 
occurs at the ends of the parameters, where it corresponds to a real infrared divergence. 
A regulator for the latter, for example, a small mass for the massless particle, also offers 
a regulator for a — > oo. 

3.5 Rotation in spatial loop momentum 

Finally, we discuss another prescription that seems to have not been applied in ordinary 
field theory calculations. As it leads to the same result as the above prescriptions for loop 
integrals in ordinary theory, it might also offer a possible means in NC theory. We first 
finish the n — 1 dimensional angular integral, 

Q roc roo k n 2 

iB^(m 2 -ie) = — — ^- / dk / dk—, —7 (29) 

5 y ' (2tt)™ 7-oc Jo (ifeg - k 2 - m 2 + it) 8 v 1 

with k = \k\. The integrand of the k integral, considered as a function on the complex 
k plane, has two poles of order 5 in the first and third quadrant. Thus, we may rotate 
clockwise the integration path from the real positive axis to the negative imaginary axis, 
k — > e~ l %k, 

iBUm 2 -ie) = jfe=l e -*f I'" dk r dk — 2 - 9 fc "" 2 9 , . . (30) 
dK ' (27r)™ J-oo Jo (fcg + k 2 - m 2 + te) 5 
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Now write it back to an n dimensional Euclidean integral, 



iBg(m 2 — ie) 



r d n k E 



1 



= e-^K"- 1 ) 



J (2k)- (kl 

^ i r 



(27r) n 2 Jo 



dx 



m 2 + ie) 5 



X 2 



n-2 



(31) 



(x — m 2 + ie) 5 



It can be worked out to reproduce the result given by the other three prescriptions. 

Note the k factor in the numerator in eq. (29). The rotation looks awkward when n 
is large enough to deny the justification of ignoring the contribution from the quarter- 
circle while Wick rotation in k still works. In such a case, we are not really considering 
a rotation. Instead, we are considering whether it is still possible to mimic the result 

— * 

obtained in the other prescriptions by integrating directly over the imaginary axis of \k\. 
Using the conception of rotation helps us find which half of the imaginary axis would be 
appropriate for this. 

3.6 Summary 

It is interesting that the above three manipulations are equivalent in ordinary theory to 
the defining one implied in Feynman propagator as far as the loop integral calculation is 
concerned. This feature is due to the simple rational structure of the integrand which has 
only the pole singularity and decays as a power for large loop momentum in all regions 
except the one close to the light-like domain. Beyond the ordinary theory, there is no 
guarantee that the three manipulations are still justified and should be equivalent to the 
defining prescription even if they look mathematically meaningful when applied by brute 
force. But it is clear from the above discussion that, as long as the manifestly covariant 
Feynman propagator appears as an ingredient of the theory, for consistency reasons we 
should stick to its fundamental defining prescription when evaluating loop integrals. This 
is the spirit that we will follow in the next section. Whether the derived prescriptions are 
still justified or equivalent arises as a curious subject that we will treat later on. 

4 Minkowski spacetime: noncommutative theory 

The loop integral of an NP diagram becomes more involved in Minkowski spacetime due 
to the indefinite metric which makes the calculations in Euclidean spacetime not directly 
applicable. In this section, we describe our calculations based on the defining prescription 
of Feynman propagator and following the lessons we learned in Euclidean spacetime. The 
potential problem discussed in subsection 3.1 is avoided as in ordinary theory. 
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4.1 Kinematic analysis 

As in Euclidean spacetime, the loop integral of an NP diagram in n dimensional Minkowski 
spacetime can be reduced to the following form, 

™/ 9 9 x f d n k cos(k-p) . n . . 

^?(P 2 , - lC ) = y ( 27r )n( fc 2_ m 2 + ^ » n ^ 2 ' mte g er ^ > 1 ? e = 0+ 

(32) 

But now p 2 can take positive (time-like), zero (light-like) and negative (space-like) values. 
For instance, for pure space-space noncommutativity, p 2 is negative semi-definite; if time 
is involved, p 2 can take either sign. 

As we shall see below, the convergence property of the integral Bg differs according 
to the sign of p 2 . Thus, to prepare for the later calculation, we first rewrite the integral 
in a standard form. A change of integration variables, a Lorentz transformation L on k 
in particular, does not change the value of the integral. This leaves the denominator and 
integration measure invariant and can be absorbed into the inverse transformation of p in 
the phase factor. Note that this may cause a problem for a nonconvergent integral when 
its regularized version is not Lorentz invariant. This will be discussed later. Therefore, 
by choosing L properly, a p with p 2 > can be transformed into a new one with p — 
sign(p )v / P ? 7^ 0, p = 0, and a p with p 2 < can be transformed into one with p = 
0, p 7^ 0, |p| = \J — p 2 . For p 2 = 0, we can make its components as small as we want 
but cannot reach the zero limit with a regular Lorentz transformation, while leaving the 
value of the integral unchanged. This is the UV/IR mixing problem when the integral Bg 
without the phase does not converge. We will not consider this last case. We emphasize 
that for the above standardization we need not assume that p is a Lorentz vector. If it 
is, p 2 is the same in all frames. If it is not, it is not a problem either: all components 
of p must be specified in a given frame before integration, and then according to this 
specified p we choose our L to make p standard. The only difference for us is that a single 
calculation holds for all frames in the former case, while we have to do the calculation 
frame by frame in the latter case according to the uniquely specified p. 

As in the case of Euclidean space, Bg and Bg^ are related by a derivative in m 2 when 
both converge, or when one or both of them are regularized in a way that does not depend 
on m 2 . We will thus concentrate on £?" and discuss briefly the results for Bg. 

4.2 Space-like case: p 2 < 

The standardized integral is 
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where only the (n—1) dimensional spatial components of k are involved in the phase. Note 
that it is not possible at this stage to do integration by parts in the variable k ■ p as we did 
in method 1 on Euclidean spacetime. This is due to the indefinite metric which cannot 
guarantee vanishing of the surface terms. This is precisely the same problem we discussed 
in subsection 3.1. According to the defining prescription of Feynman propagator, the ko 
integral should be finished first using residue theorem for the prescribed contour closed 
in the upper or lower half complex plane of k , 

„ . 9 9 If d n ~ l k cos(k-p) 

iWW-) = -J^^JL. ,34) 



We emphasize that the integration order matters here. With the ko contour integral 
done first, it is understood that it is finished with respect to a fixed and finite \k\ so that 
one of the two poles stays always inside the contour. Without enforcing this, the ko and 
k subintegrals would be interwound since there would be a domain where both poles stay 
outside the contour when, roughly speaking, \k\ is larger than the radius of the contour. 
It is also in this limited sense of the integral that we can consider Wick rotation (see 
subsection 5.1) and argue meaningfully that the contribution from the semi-circle of large 
radius can be ignored when the integrand drops fast enough for large \ko\ but finite \k\. 

Now the Euclidean subtlety just mentioned sets in, as it is mathematically possible 
to do integration by parts in k ■ p and ignore the surface terms, followed by integration 
over perpendicular components, or vice versa. The two manipulations yield the best and 
the worst convergence respectively. As in the Euclidean case, we have to appeal to the 
physics context of field theory where the integral appears, to argue against manipulations 
with reference to the direction specified by the p. This is the only subtlety in calculations 
using the defining prescription and it seems unavoidable. 

The angular integration yields, 

3 — n 

B[ l (p 2 ,m 2 -le) = --(2tt) ' / dk k n ; , - 

I Jo vfc + m 2 — te 

Using eq. (8), we see that is not well-defined in the UV for n > 4. Note the difference 
from the Euclidean case concerning convergence at n — 4. There, the angular averaging 
is over n dimensions, while here it is over n — 1 dimensions since ko has been first finished. 
Generally, Bg converges for n < 45, oscillates in the UV at n = 45 and diverges oscilla- 
torily for n > 45. As naively expected, the NC phase in the case p 2 < indeed improves 
the convergence property of the loop integral compared to its commutative counterpart. 
To regularize the UV behaviour, we introduce a spherical cut-off on the spatial loop 
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momentum, \k\ < A. Rescaling k\p\ = x and denoting k = A|p|, we have, 

X~2~ Jn-3 (X) 



I i_„ /-ft X 2 J„_3{X] 

I^(p 2 ,m 2 -* e ) = --(27t)^|p1 2 ""/ rfx 2 2 , « 

^ Jo y x 1 — a z — ie 



2 2 2 

— m p 



(36) 



For clarity, we distinguish the two cases of m 2 > (i.e., a 2 < 0) and m 2 < (i.e., a 2 > 0). 



The case m 2 = can be covered by either of them in the limit 



m 



0^ since has no 



IR divergence for n > 2. From now on we specialize to four dimensions though there is 
no difficulty to work out new, more divergent terms in higher dimensions. The results for 
lower dimensions will be discussed in a separate subsection. 

4.2.1 m 2 > 

It is safe to ignore — ie since the argument of the square root is positive definite. The 
non- convergent term can be separated by integration by parts, 



i27T) 2 \p\ 2 Bf(p 2 ,m 2 -ie) 



xsmx 
dx—j= 

^jx 2 +\a 2 \ 

KCOSK 2 f K COSX 

= + \a \ ax s- 

yW|a 2 | -A) (x 2 + |a 2 |)l 



(37) 



The second term is now well-defined in the limit A — > oo. Keeping terms not vanishing 
in the limit, we have, 

1 



(27r) 2 Bf(p 2 ,m 2 -ie) = 



cos 



(38) 



4.2.2 m 2 < 

Since a 2 > 0, the argument of the square root can become negative for small enough x so 
that we must keep — ie in that region. Note that the integral is well-defined at x — \/~a?. 
For the sake of safety, we split the integral into two parts, 



-(2ir) 2 \p\ 2 Bf(p 2 ,m 2 -ie 



/•vr 

* / 

JO 



x sin x 
ax f= — == + 



x. 



x sm x 



_ dx 

a? \Jx 2 - a 2 



(39) 



7T 



The first integral equals i—Va 2 Ji yva 2 ). The second one is not well-defined in the UV 
and can be worked out using the following trick, 

xsinx r K x 



f 



dx 



\/a? \J x 2 — a 2 



\/a? ^ x 2 — a 2 
K 



d 



cos x — cos V cr 



\/k 2 — a 2 



cos k — cos V a 



dx 



cos x — cos v cr 
(x 2 — a 2 )i 



(40) 
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The remaining integral now converges and equals in the limit A — > oo, 

— — Va?Ni (^faPj — cos Va?, 

where N v {z) is the Bessel function of the second kind of order v. Collecting terms and 
keeping only terms that survive the limit, we have 

1 



-{27r) 2 Bf{p 2 ,m 2 - ie) 



P 2 



. , cos k — i—\a l 
2 1 2 



Ji (v 7 ^) + iN Y (v 7 ^) 



cos 



K - (-iV^)^ (- 



(41) 



where a relation amongst Bessel functions has been used. 
The results for space-like p are summarized compactly as 



i2n) 2 Bf(p 2 ,m 2 - ie) 



1 



cos [AJ-p 2 ) - zqK^zo) 



\2v) 2 B\{p 2 ,m 2 -ie) = --K (z ), 



(42) 



with z = ^ —p 2 (m 2 — it). Note that the first term in Bf oscillates in the UV though it 
is bounded. It is missed in the conventional calculations. B\ and the second term in Bf 
coincide with the conventional results. 

4.3 Time-like case: p 2 > 

The standardized integral becomes in this case, 

d n k e ikopo 



iB™(p 2 ,m 2 -ie) = J 



— , p = VP > 



(43) 



(2?r)" k 2 - m 2 + ie' 

Note that we have used the reflection symmetry to replace cos(/coPo) by the phase in the 
sense of principal value. The same symmetry in ko allows us to take either sign of po- This 
affects the choice of contour for the k integral accordingly: for each choice of sign(p ), 
only one of the two defining contours makes sense; but either choice leads to the same 
answer of course. For the above choice, we close the contour in the upper half complex 
plane of k . Picking up the residue of the pole in the second quadrant and finishing the 
trivial angular integration yields, 

-ipoVk 2 +m 2 — ie 



B™(p 2 , m 2 — ie) = 



1 n 



n-1 



2(2tt) 



n-1 



dk k 



n-2 



\Jk 2 + m? 



(44) 



ie 



Taking into account the oscillating phase, the above integral is not well-defined for n > 3. 
Note that there is no similar Euclidean subtlety as occurring in the case p 2 < 0. The 
above integral is uniquely specified. One might argue to do integration by parts and ignore 



17 



surface terms. This would make any theory free of divergences as it simply amounts to 
ignoring any divergences. For the same |p 2 |, the integral in the time-like case is less 
convergent than in the space-like case. This arises from the asymmetry in the number of 
separate dimensions in time and space: while the oscillating behaviour in the radius k is 
similar, the angular oscillation occurring in the space-like is lacking in the time-like case. 
This asymmetry disappears in two dimensions. As above, we introduce a cut-off on the 
spatial loop momentum k < A to regularize the UV behaviour. Rescaling kp = x and 
denoting k = A^/p 1 , we have, 



B[ l (p 2 , m 2 — it) = 



1 n 



n-l 



2-n 



dx 



x'- 



-2 „— W x 2 +a 2 —ie 



2 J2 



V X 1 + a 2 — 
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m p 



(45) 



2 (27T)"- 1 

Below we work out the integral in four dimensions according to the sign of m 2 . For a gen- 
eral discussion concerning convergence of Bg, see subsection 4.3.3. For lower dimensional 
results, see subsection 4.4. 

4.3.1 m 2 > 

In this case a 2 > 0, the integral reduces in four dimensions to 



-(2n) 2 p 2 Bf(p 2 ,m 2 - ie) = T 

Jo 



dx 



= ine 



\Jx 2 + a 2 

-iVK 2 +a 2 _ ■ r 

Jo 



dx e ~ tV7IT ^ 

The remaining integral is further reduced in UV divergence as follows, 



(46) 



dx e 



= %- 



dt 



.VK 2 + a 2 



Vt 2 - a 2 



K 



^-iVn 2 +a 2 _ ^-iVa 2 



+ ia 



Vn 2 +a 2 e -it _ e -iVa- 

21 dt 



(t 2 - a 2 



(47) 



The last term is now well-defined and equals in the limit A — > oo, 
Thus, we find, 



in + -a A + 1 ) e lK - iV a 2 



i2n) 2 Bf(p 2 ,m 2 -ie) = - 

f 

i2n) 2 Bl(p 2 ,m 2 -te) = -[e' iK - K (iVtf) 



K 1 (iVr^ 2 ) 



(48) 
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4.3.2 m 2 < 

As we did with the case p 2 < 0, m 2 < 0, we split the integral into two parts: 



-{2n) 2 p 2 Bf{p 2 ,m 2 - ie) 
The first term equals 



L 



dx 



+ 



cr — ar 



dx — , (49) 



x^ — cr 



.7T 



h 



Li Via 



where /^(-z) is the other modified Bessel function and L u (z) is the modified Struve function. 
The second integral is reduced by integration by parts until we end up with a convergent 
one in the limit k — > oo: 



f 



. dx 



x 2 e -iy/x*-\a*\ 



= a 



ar — a z 



i^ 2 \V^Ti + 



dt V^Tle- u V\ a2 \ 
t 



ity/\a?\ 



k 2 l 



(50) 



Jo 



where the last integral is well-defined in the UV and equals 

71 



-\\a 2 \K x 



i-\l\a- 



L-A\\a?\\-h 



Collecting terms, keeping only those that survive the limit and using L_i(z) — Li(z) = — 

TV 

to cancel one of the surface terms, we have, 



-(2n) 2 Bf(p 2 ,m 2 -te) = 
-(2n) 2 Bl(p 2 ,m 2 - ie) = 



IK 



-K n 



o V a 



(51) 



The results for time-like p are summarized as 

iA 1 



-(2n) 2 Bf(p 2 ,m 2 - ie) = 
-(2ir) 2 Bl(p 2 ,m 2 - ie) = 



'P 



1 



ZqK^Zq), 



e- iA ^-K (z ) 



(52) 



with zq = y—p 2 (m 2 — ie) again. They are indeed very different from those in the con- 
ventional approaches: B\ oscillatorily diverges in the UV while B\ oscillates in the UV, 
though the K function terms coincide. It is interesting to note that these new terms are 
complex, independently of the sign of m 2 . 
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4.3.3 Improving or worsening convergence? 



We saw in the above that the NC phase in the case p 2 > softens an otherwise quadrati- 
cally divergent integral to an oscillating linearly divergent one and an otherwise logarith- 
mically divergent integral to an oscillating but bounded one. However this improvement 
of convergence is not a general tendency, but a special feature of the integral To see 
this, we consider the integral Bg for general integers n and 5 > 1 in the case p 2 > 0. 
According to the defining prescription, we finish first the k contour integral by residue 
theorem, 

1 



B$(p 2 ,m 2 -ie) = 



(27T)™- 1 T(5) Jo 



x 



dk k 



n-2 



d: 



8-1 



exp [ik{ 



(53) 



dk 5 1 (k - Vk 2 + m 2 - ie) s \ fc0= „ V fc . +m ._ t£ 

Now the derivative can act on the denominator and the phase. When it acts on 
the phase, it does not lower the divergence degree in k so that we lose in convergence 
compared to the ordinary theory. This is a consequence of the interplay that we referred 
to at the beginning of section 3 between the complexation of the phase and the propagator 
on the ko plane, which makes Wick rotation impossible and causes the result from the 
contour integration less convergent. It seems inevitable for the correct manipulation of 
propagators to avoid the apparent problem in subsection 3.1. The fate of convergence 
thus rests on the term where all derivatives act on the phase, which is 



5-1 



(27T)"- 1 (-2)*r(5) Jo 



f 

Jo 



dk k 



n-2 



exp 



isjp 2 {k 2 + 



77T 



16, 



+ 



(54) 



[yk 2 + m 2 — iej 

where the dots stand for less divergent terms. Using a cut-off A, the leading divergence is 



f>n 

n s — 



(27T)™- 1 (-2YT(5) v 7 ^ 



A 



+ 



(55) 



Thus the integral converges for n < 5 + 2, oscillates at n = 5 + 2 and diverges oscillatorily 
for n > 5 + 2. For comparison, the ordinary integral without a phase converges for 
n < 25, diverges logarithmically at n = 25 and as a power A n ~ 2S for n > 25, with the 
same UV cut-off. Thus, only at 5 = 1, the NC phase improves the convergence property 
of the integral by reducing the divergence degree by one. At 5 = 2, the divergence in 
ordinary theory is modulated by an oscillating phase, except at n = 4 where a logarithmic 
divergence softens to a bounded oscillating factor. For 5 > 3, convergence can be worsened 
in high enough dimensions; i.e., a convergent integral in ordinary theory (when n < 25) 
is rendered divergent (n > 5 + 2) or oscillating (n = 5 + 2) by the phase. For 5 > 1, Bg 
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contains 5 terms that successively differ by one in divergence degree. The terms whose 
convergence is improved over the commutative case are singular at p 2 = 0, which is the 
UV/IR mixing problem, while the terms whose convergence is worsened or just modulated 
by a phase are smooth in the limit p 2 — > 0. 

4.4 Results in two and three dimensions 

For completeness, we include in this subsection the results in lower dimensions. The two 
dimensional case is interesting because time and space are more symmetric than in higher 
dimensions. Here, we have 6^ = 9e^ u with e M1 , being the Lorentz invariant antisymmetric 
constant tensor, so that p 2 = —9 2 p 2 . The cases of p 2 > and p 2 < are then related by a 
sign flip of the external momentum quadrature and there should be no essential difference 
between the two. This indeed coincides with our previous results that all NP integrals 
are convergent in two dimensions. Actually, any of the four prescriptions may be used 
to evaluate the integral as long as it makes sense mathematically. In what follows, we 
establish the above mentioned relation more directly. 
It is sufficient for us to consider Bf, 



= / J^f k 2 - rn? + ie * ^ 



cos(A;ov / /?) for p 2 > 
cos(/c lv / — p 2 ) for p 2 < 



(56) 



We relate it to Bf(—p 2 , —m 2 —ie). Complex conjugation and interchanging the integration 
variables for the case p 2 > yields, 



d 2 k cos(k lx /-(-p 2 

—iB{~{p",m" — ie) = — J - 

namely, 



-zB 2 *(o 2 m 2 -te) - I d k —{^-y-P)J (57) 
IB, [p,m it) - J ^ fc2 _ ( _ m2)+ . e , 



Bf*(p 2 ,m 2 - ie) = Bf(-p 2 , -m 2 - ie) (58) 

which now holds for any sign of p 2 . Thus, only the relative sign of p 2 and m 2 is relevant. 
The explicit result is 

-(27t)B 2 (p 2 ,m 2 - ie) = K (z ) (59) 

Thus B\ is real for p 2 m 2 < and complex for p 2 m 2 > 0. 
The results in three dimensions are 

v£ D 3, 2 2 Vm 2 - ie T ^ , x , [ _ forp 2 <0 



-(2n)-*Bt(p\m 2 -ie) = — J,(, ) + , & for p 2 > ^ 
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4.5 Further discussions 

Before we started off to evaluate the integral, we first made the integral standardized. 
This is all right if it is well-defined, as it amounts to a change of integration variables. 
For integrals that are not well-defined however, this change of variables may clash with 
the regularization used to define them. This is indeed the case with the cut-off methods 
as we discussed at the end of subsection 2.2 for the Euclidean case. Since a suitable reg- 
ularization should maintain Lorentz symmetry (suppose for simplicity that p is a Lorentz 
vector), all reference to individual components of p must be removed before the cut-off is 
imposed. This is the other reason that we first standardized the integrals, in addition to 
the one explained in subsection 4.1. But differently from the Euclidean space, we cannot 
regularize the UV behaviour by simply imposing \k 2 \ < A 2 due to the apparent problem 
exposed in section 3.1. In ordinary theory this is circumvented by making a Wick rotation 
beforehand. However, as we explained previously and will detail in the next section, there 
is generally no guarantee for the justification of Wick rotation in NC theory. Furthermore, 
in the above approach based on the defining prescription, we never put a cut-off on k as 
this is demanded by self-consistency of the Feynman propagator; instead, if necessary, the 
cut-off is always imposed on spatial components after finishing the ko contour integration. 
We stress that the new terms in eqs. (42, 52) compared to the conventional results are 
not regular in the UV and do not disappear in the limit of A — > oo. 

5 Ambiguities on noncommutative Minkowski space- 
time 

In the preceding section we evaluated the NP integrals according to the defining prescrip- 
tion built in Feynman propagator. This is the fundamental prescription that is justified 
by consistency whenever Feynman propagator arises as an ingredient of any theory. In 
that evaluation, only one subtlety remains for p 2 < 0, which is essentially the same one as 
in Euclidean spacetime exposed in section 2. As it cannot be excluded mathematically, it 
has to be fixed by physical motivations. The results are thus certain up to this subtlety. 
Note that the latter appears both for space-space and time-space noncommutativity. 

In this section, we employ the other three prescriptions to calculate the NP integrals 
when this seems to make sense mathematically. Our purpose is to show that equivalent 
manipulations for ordinary theory are not necessarily equivalent any more in NC theory. 
Ambiguities arise in the sense that these three prescriptions generally lead to different 
results concerning convergence than the defining one. In the course of doing so, it will 
also become clear how the difference occurs between the calculation in the preceding 
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section and the conventional approaches. 
5.1 Wick rotation 

When p 2 < 0, k is not involved in the NC phase in the integral B™. It thus makes 
sense mathematically to integrate over the imaginary axis of k , which amounts to Wick 
rotation in fc , 



The above integral taken in its face is not uniquely specified as different sequential ma- 
nipulations can result in different answers concerning convergence. We will not go into 
the details but listing the results. (1) If we finish k n first, the result coincides with that 
in the preceding section, including the same Euclidean subtlety. This is just as expected, 
since we know that ko before rotation, corresponding to k n after rotation, should be fin- 
ished first. The coincidence is guaranteed by the analytic properties. The following three 
manipulations are similar to the Euclidean case. (2) If we first do n dimensional angular 
averaging treating p as a vector in n dimensions, the integral converges for n < 5. This 
looks not quite natural. Why should we still do an n dimensional spherical averaging 

— # 

now that p has no n-th component? (3) If the component k ■ p is finished first, the inte- 
gral converges for any n. (4) Since k\ > 0, the denominator may be exponentiated by 
Schwinger parametrization. Even if m 2 < 0, this is still all right since the k E domain 
with the denominator negative is finite and thus the exponentiation cannot modify the 
UV property in k E . This is a popular approach in the literature. However, as explained 
in subsection 2.1, for the parametrization to go through, we have to assume that the 
integral converges. It is thus not surprising that its outcome converges for any n in this 
manipulation. 

For p 2 > 0, there is no way to invoke Wick rotation since the NC phase blows up on 
either the positive or negative imaginary axis, depending on the choice of the phase sign. 
In Ref. [3], this is circumvented by simultaneous analytic continuation of the external 
momenta p and 6^ v (and thus m 2 and p) together with k , so that a phase on the real 
axis remains to be a phase on the imaginary axis. [Note that 'simultaneous' continuation 
is vital to avoid blowing-up exponentials. There is no step in between, since it does not 
help to continue p and 6^ v before or after continuing k . Such a continuation is very 
unusual indeed. It was also used in Ref. [3] to argue away a potential UV divergence 
in NP loop integrals.] After finishing the integration, the external momenta and M „ are 
continued back to Minkowski spacetime. However, it is precisely this continuation forth 




(61) 



where k^ n ^ reminds us that the rotated n-th component is not involved in the phase. 



23 



and back in p and 6^ that has to be justified, instead of being assumed. We know from 
the above discussion that convergence properties on NC Euclidean spacetime are very 
different from those on NC Minkowski spacetime for p 2 > due to their big difference 
in analytic properties: in the former, everything else except the phase is positive-definite 
and thus the integral can be finished in principle without continuation from the real axis 
to the complex plane; in the latter, analytic continuation is necessary due to the apparent 
problem associated with the indefinite metric, which however is made nontrivial by the 
NC phase. 

5.2 e parametrization 

This manipulation is more tricky and looks harder to exclude than the other two. It starts 
with the justified Schwinger parametrization based on the infinitesimal e > 0: 

i 5+1 T{5)B%{p 2 ,m 2 -ie) = f (£y^ j™ da a *-i e -«H(* 2 — 2 )] , (62) 

with integer 5 > 1. As we will show below, this only postpones the problem instead of 
solving it. To proceed further, we have to interchange the integration order of k and a 
so that we can finish the loop integral. As we stressed repeatedly, this is legitimate only 
when the original integral is well-defined. To see that the problem will come back, let us 
tolerate this and press on. Complete the k quadrature and make a shift k k — p/{2a). 
The k integral is finished as in the ordinary case, 

J o 

Since the a -1 term in the exponent is purely imaginary, the above integral is not well 
behaved at a — > for n > 2(5 + 1) and hence cannot be directly evaluated in terms of 
Bessel functions. This is the standard UV divergence though seemingly improved, but 
it implies already that the NP integral Bg is not always UV finite. Furthermore, as in 
the ordinary case, e is a routing indicator whose absolute magnitude should be irrelevant. 
Namely, to justify the parametrization posteriorly, dependence of convergence on e at 
a — > oo should be relaxed. We recall that this is the region which causes an apparent 
problem due to the indefinite metric in the ordinary case. Nevertheless, for the most 
frequently studied diagram in the literature with two propagators in four dimensions, i.e., 
.Bf , the above seems to give a well-defined result even if we set e = 0, 

poo /}2 

(47r) 2 r(2)5 2 4 (p 2 ,m 2 ) = / da a" 1 exp -iam 2 - i^- 

2K (sign(p 2 ) iy/pW) for p 2 m 2 > ( 64 ) 



. / 2 ■ \ .9 
—ta(m — le) — %-— 

4a 



(63) 
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2K (yV™ 2 !) for P 2 m 2 < 
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which may be summarized concisely as the conventional result, 



(47r) 2 r(2)£ 2 V,m 2 ) = 2*o 



> 2 {m 2 — ie) 



(65) 



However, the above manipulation is rather questionable. We started with a positive e 
whose non-vanishing is crucial for us to apply the Schwinger parametrization and whose 
role waits to be relaxed. But finally we found it is redundant for the convergence of the 
parameter integral. Although it appears in the final compact result, this re-introduction 
is only for the purpose of bookkeeping. To see the point more definitely, let us apply the 
same practice of setting e = to the integral _B|, 



safest prescription that B\ converges for either sign of p 2 . Thus, the role played by a 
nonvanishing e cannot simply be ignored. The only way left for the feasibility of the 
method is to relax the dependence of convergence on e by rotating a from the positive 
real axis to the positive or negative imaginary axis. 

We showed in subsection 3.4 that the relaxation is always possible in the ordinary case, 
which thus confirms the equivalence of the parametrization to the defining prescription. 
In the NC case however, this is not always possible. For p 2 m 2 > (i.e., below threshold 
for p 2 > or above threshold for p 2 < 0), the a and a~ l terms have an opposite sign in 
the real parts for complex a. It is thus not possible to rotate it to either its positive or 
negative imaginary axis, to avoid exponential blowing-up. Note that the power factor in 
eq. (63) does not help on this matter. This is strong enough to deny the feasibility of 
the approach, though it seems to work for p 2 m 2 < 0, because, for each sign of p 2 both 
kinematic configurations of m 2 > and m 2 < are permitted and cannot be covered 
simultaneously. 

5.3 Rotation in \k\ 

For p 2 < 0, analytic continuation in \k\ is not possible: since we are only allowed to rotate 
it in the fourth quadrant due to its pole in the first quadrant, this rotation is allowed only 
for the e~ l '^v^~^ term but not for the e +l ^^~^ term, both of which appear after angular 
integration. For p 2 > 0, the rotation is possible mathematically, up to the explanation in 
subsection 3.5. Consider again B™, which becomes after angular integration, 




(66) 



which obviously diverges at a 



oo. But we know from the last section using the 



iB™(p 2 , m 2 — ie) 




dk cos(k p ) 



k n ~ 2 dk 



(67) 
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o kn — k 2 — m? + ie 
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We rotate k in the fourth quadrant, k — > e 1 2 k, to move to the Euclidean spacetime; then 
we put back the angular integration, 

Dn/ 2 2 • \ -j^TT f dnk E cos(A; a//?) ( . 

^,m'-«) = e . J— kl _ m2 + u (68) 

Note the opposite sign in the denominator compared to Wick rotation. 
Rescaling kE\fp 1 = x and relabelling x, the integral becomes, 

WW" 1 -') = + u < 69 > 

As far as the UV property is concerned, the above integral is the same as that we treated 
in Euclidean spacetime. This is obvious for m 2 < 0. When m 2 > 0, the denominator has 
a sign flip at a finite value of x, which does not modify the UV property and is taken 
of by +ie. The integral can be similarly worked out with same subtleties as appearing 
there, which we do not repeat. For comparison, we recall that there is no subtlety for 
p 2 > in the evaluation using the defining prescription. All this shows clearly that the 
prescription of rotation in \k\ is no more equivalent to the defining one in NC theory. The 
weakest point of rotation in \k\ is perhaps its lack of physical motivation compared to 
Wick rotation, which in ordinary theory relates Minkowski to Euclidean spacetime. 

6 Application to real Lp A theory 

To illustrate physical implications of the above calculations, we consider briefly the real 
(/? 4 theory in four dimensions, 

C = —d^ip -k d^ip m 2 (p-k(p (p * ip * ip ★ (p (70) 

2 2 4 

The Feynman rule for the vertex is 

-i2\[c 12 c 34 + c 23 ci4 + C31C24] (71) 

with Cij = cos(pi A pj), where pi is the incoming momentum of the i-th particle, and 
P A q = OuutpPq". 

The NP part of the one- loop self-energy, defined as ix diagram, is 

S NP (p 2 ,m 2 -ie) = -\B*(p 2 ,m 2 -ie), (72) 

where p is the external momentum. This diagram is usually not checked probably because 
it is too simple. But according to eqs. (42,52), it is not. For p 2 > 0, which is possible 
when 9 j 7^ 0, the standard term proportional to — zqKi(zq), with z = yj ' —p 2 {m 2 — ie) = 
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iy/p 2 m 2 , violates unitarity, though its planar counterpart is empty in unitarity. The new 
terms are oscillatorily divergent in the UV. Since they are complex, they violate unitarity 
too, and furthermore the violation is even not finite. For p 2 < 0, the diagram has a 
new term which is bounded but oscillates in the UV. Note that the new terms cannot be 
absorbed by renormalization of the bilinear terms in C. 

Consider now the one-loop vertex with incoming external momenta pj. There are three 
diagrams related by crossing symmetry. It is sufficient to consider the contribution from 
the diagram where particles 1 and 2 merge into 3 and 4. After some tedious algebra, the 

(12) 

NP part of the vertex, T NP , defined as ix diagram, is reduced to 





(pI 


A 12 - ie) 


+ Bt 


(pI 


A12 - ie) 




(pI 


A 12 - ie) 


+ Bt 


{pI 


A 12 - iej 



1 

-^2 r NP (Pi,Pj) = J dx {2c 12 c u B2 (pu, A 12 - ie 

+C12 cos((2x - l)p 3 Ap 4 ) 
+c 34 cos((2a; - l)p 1 Ap 2 ) 

+- cos((2:r - l)( Pl A p 2 - p 3 A p 4 ))B* (p 2 13 , A l2 - it) 
f 

+- cos((2a; - l)( Pl A p 2 + p 3 A p A ))B% (p 2 u , A 12 - ie) } 

(73) 

where A X2 = m 2 — p\ 2 x{l — x), p^ a = 6^ v p v a with a = 1, 2, 3, 4, 12, 13, 14, and p X2 = pi + p 2 , 
etc. When all p 2 < 0, as is the case with pure space-space noncommutativity, there is 
nothing new compared to the conventional results. When some p 2 > however, there 
will be new contributions to the vertex. We will not attempt here a detailed kinematic 
analysis to exhaust all possibilities but giving a few examples. For instance, the separate 
new term from p\ 2 > 0, p\ > and pf 3 > is 



-(2n)- 2 c 12 c M e~ lA V&, 

I, sin(p 3 Ap 4 ) „-iA^ftf 
--(2tt) C12 e Vi, / 74 n 

2 p 3 A P4 v ' 

1 (277)~ 2 Sil1 ^ 1 ~ P3 AP4 " > c - i ^.^ 

4 Pi A p 2 - p 3 A p 4 

Each of these terms oscillates in the UV, is complex thus harming unitarity, and smooth 
in the commutative limit. 



7 Conclusions 

We have studied the loop integrals appearing in one loop NP diagrams of NC field theory 
and obtained results that are very different from those in the literature. We found that a 
mathematical subtlety occurs already on Euclidean spacetime. We appealed to a physical 
motivation in the context of quantum field theory to fix it up. With this manipulation, 
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NP integrals converge in four and lower dimensions but can diverge or oscillate in the 
UV in higher dimensions. We employed a simple cut-off method to illustrate the new UV 
non-regular terms that were missed in the conventional approaches. 

The situation becomes more involved on Minkowski spacetime due to an interplay 
between the indefinite metric and complexation of the NC phase. For consistency rea- 
sons that positive (negative) frequencies must propagate forward (backward) in time, the 
manifestly covariant Feynman propagator has to be defined as a contour integral on the 
complex plane. This avoids an apparent problem due to the indefinite metric in ordinary 
loop integrals. We emphasized that it is thus mandatory to evaluate NP integrals of NC 
theory in this defining prescription where the zero component of the loop momentum is 
first worked out by the contour integration. It has to be examined as a separate issue 
whether other prescriptions that are equivalent in ordinary theory yield the same result as 
the defining one. Indeed, an NC phase defined on the real axis ceases to be a phase on the 
complex plane, but can blow up in the upper or lower half plane when the NC external 
momentum is time-like, for instance. This directly blocks up the usage of some of the 
other prescriptions in the NC case. The remaining prescriptions, even if they may seem 
to work mathematically, are either not physically well-motivated (e.g., rotation in \k\ in 
the case of a time-like external NC momentum), or cannot go through for all possible 
kinematic configurations (e.g., e parametrization), or can be argued away using the same 
physical motivation we appealed to on Euclidean spacetime (e.g., Wick rotation for a 
space-like NC external momentum). These discussions show that equivalent prescriptions 
in ordinary theory are not necessarily equivalent in NC theory. 

An NC phase involving a space-like NC external momentum always improves upon 
the convergence property. In four dimensions, such a phase softens an otherwise quadrat- 
ically divergent integral to a bounded and UV oscillating one, while it saves a logarithmic 
divergence to a finite result (up to the UV/IR mixing as always). In higher dimensions 
however, new UV divergent terms appear that are modulated by an oscillating factor. 
When a time-like NC external momentum is involved, the NC phase does not always im- 
prove convergence but can worsen it in high enough dimensions. This is a surprising result 
indeed and has its origin in the interplay discussed above. It lends a strong support to the 
opinion that Feynman propagator does not appear as a basic ingredient in perturbative 
NC theory when time does not commute with space [8, 9, 10] so that the naive approach 
built on it is not well-founded. In four dimensions, such a phase softens a quadratically 
divergent integral to a linearly divergent one modulated by an oscillating factor, and a 
logarithmically divergent integral to an oscillating one. We determined those new terms 
using the cut-off regularization in the defining prescription. As we demonstrated in the 
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real <p 4 theory, they endanger renormalizability already at one loop level, and become 
complex in the case of a time-like NC external momentum and thus harm unitarity as 
well. 
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